Fermionic Casimir densities in anti-de Sitter spacetime 

E. Elizalde^: S. D. Odintsov^'^t A. A. Saharian^^ 

^Instituto de Ciencias del Espacio (CSIC) 
and Institut d'Estudis Espacials de Catalunya (lEEC/CSIC) 
Campus UAB, Facultat de Ciencies, Torre C5-Parell-2a planta, 
08193 Bellaterra (Barcelona) Spain 

"^Institucio Catalana de Recerca i Estudis Avangats (ICREA) 

^Department of Physics, Yerevan State University, 
1 Alex Manoogian Street, 0025 Yerevan, Armenia 

February 13, 2013 



Abstract 

The fermionic condensate and vacuum expectation value of the energy-momentum tensor, for 
a massive fermionic field on the background of anti-de Sitter spacetime, in the geometry of two 
parallel boundaries with bag boundary conditions, are investigated. Vacuum expectation values, 
expressed as series involving the eigenvalues of the radial quantum number, are neatly decomposed 
into boundary-free, single-boundary-induced, and second-boundary-induced parts, with the help of 
the generalized Abel-Plana summation formula. In this way, the renormalization procedure is very 
conveniently reduced to the one corresponding to boundary-free AdS spacetime. The boundary- 
induced contributions to the fermionic condensate and to the vacuum expectation value of the 
energy density are proven to be everywhere negative. The vacuum expectation values are expo- 
nentially suppressed at distances from the boundaries much larger than the curvature radius of 
the AdS space. Near the boundaries, effects related with the curvature of the background space- 
time are shown to be subdominant and, to leading order, all known results for boundaries in the 
Minkowski bulk are recovered. Zeta function techniques are successfully used for the evaluation 
of the total vacuum energy in the region between the boundaries. It is proven that the resulting 
interaction forces between them are attractive and that, for large separations, they also decay expo- 
nentially. Finally, our results are extended and explicitly translated to fermionic Casimir densities 
in braneworld scenarios of Randall-Sundrum type. 



1 Introduction 

Anti-de Sitter (AdS) spacetime is among the most popular background geometries in quantum field 
theory. Much of the earlier interest in this geometry was motivated by questions of principal nature, 
mainly related with the quantization of fields on curved backgrounds. The AdS spacetime has maximal 
symmetry and, because of this, numerous physical problems can be exactly solved on this background. 
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The presence of both regular and irregular modes and the possibility of getting an interesting causal 
structure lead to several new and remarkable phenomena. Further interest in this subject arose from 
the discovery that the AdS spacetime generically arises as a ground state in extended supergravity 
and string theories, what is again potentially most important. 

In recent developments of the topic, the AdS geometry is an arena for two classes of models. 
The first is the AdS/CFT correspondence (for a review see [I]), which represents a realization of the 
holographic principle and relates string theories or supergravity in the AdS bulk with a conformal field 
theory living on its boundary. This correspondence has many interesting consequences and provides 
a powerful tool for the investigation of physical effects in gauge theories. The second class of models 
with the AdS spacetime as background geometry is a realization of the braneworld scenario with 
large extra dimensions and provides a solution to the hierarchy problem which arises between the 
gravitational and electroweak mass scales (for reviews on braneworld gravity and cosmology see [2]). 
Here the small coupling of 4-dimensional gravity is generated by the large physical volume of extra 
dimensions. Braneworlds naturally appear in the string/M theory context and provide a novel set up 
for the discussion of phenomenological and cosmological issues related with extra dimensions. 

In the present paper, as an example of an exactly solvable physical problem in AdS spacetime, 
we will consider the Casimir effect (for reviews see [3]) for a fermionic field obeying bag boundary 
conditions on two parallel plates. The explicit dependence of the characteristics of the vacuum on 
the geometry of the background spacetime is among the most interesting topics in the investigation 
of the Casimir effect. As usual, all relevant information is encoded in the spectrum of the vacuum 
fluctuations and, not surprisingly, analytic solutions can only be found, in general, for highly symmetric 
geometries. Specifically, the Casimir effect for a massive scalar field with general curvature coupling 
parameter in the geometry of fiat and spherical boundaries on the background of de Sitter spacetime 
has been investigated recently in [1] and [5], respectively. 

Investigations of the Casimir effect in AdS spacetime have already attracted a great deal of atten- 
tion, motivated by Randall-Sundrum type braneworld scenarios [6j. In these models the background 
solution consists of two parallel fiat branes embedded in a 5-dimensional AdS bulk. The fifth coordi- 
nate is compactified on /Z"^, and the branes are on two fixed points. The fields which propagate in 
the bulk will give Casimir-type contributions to the vacuum energy and, as a result, to the vacuum 
forces acting on the branes. The Casimir effect provides in this context a natural mechanism for 
stabilizing the radion field, as required for a complete solution of the hierarchy problem. In addition, 
the Casimir energy gives a contribution to both the brane and the bulk cosmological constants and, 
hence, it has to be taken into account in any self-consistent formulation of the braneworld dynamics. 
The Casimir energy and the corresponding forces for two parallel branes in AdS spacetime have been 
evaluated in Refs. [71|8], both for scalar and fermionic fields, by using either dimensional or zeta func- 
tion regularization methods. Local Casimir densities were considered in Refs. [QtllOtlTT]. The Casimir 
effect in higher-dimensional generalizations of the AdS spacetime with compact internal spaces has 
been investigated in |12j while the Casimir energy for a massless fermionic field with generalized bag 
boundary conditions in 3-dimensional AdS spacetime was discussed in [13], for the geometry where 
one of the boundaries coincides with the AdS boundary. 

In braneworld models, two distinct types of boundary conditions arise for Dirac fermion fields, 
corresponding respectively to even and odd fields (untwisted and twisted boundary conditions) [I4] . 
In the present paper, we investigate one-loop quantum effects for a fermionic field in AdS spacetime, 
induced by two parallel boundaries with bag boundary condition. Although this condition is different 
from those appearing in braneworld scenarios, we will here show in detail how, from our formulas, the 
corresponding results for braneworlds are readily obtained. The important quantities that characterize 
the local properties of the fermionic vacuum are the fermionic condensate (FC) and the vacuum 
expectation value (VEV) of the energy-momentum tensor. For the investigation of these quantities 
we use the direct mode summation approach, which requires the knowledge of a complete set of mode 
functions for the fermionic field obeying the boundary conditions. 
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In the next section, we describe the geometry of the problem and construct the corresponding 
mode functions. By making use of them, in Sect. [3] we evaluate the FC. Applying the generalized 
Abel-Plana formula, the FC is decomposed into three parts: a boundary-free contribution, a single- 
boundary-induced one, and a second-boundary-induced one. The behavior of the VEVs in asymptotic 
regions of the parameters is discussed. In Sect. [4] we present similar considerations for the VEV of the 
energy-momentum tensor. Casimir forces acting on the boundaries and the corresponding Casimir 
energy are investigated in Sect. \5\ The generalization of the formulas obtained to the important case 
of a bulk fermionic field in the Randall-Sundrum braneworld model is discussed. Finally, our main 
results are summarized in Sect. [H 

2 Fermionic mode functions 

We consider a quantum fermionic field, in (D + l)-dimensional anti-de Sitter spacetime, AdSn+i- 
In Poincare coordinates, the corresponding line element reads 

ds"^ = gf.^dxf'dx" = e-'^y/'^riikdx'dx^ - dy^, (2.1) 

with rjik = diag(l, —1, . . . , —1), i,k = 0, . . . , D — 1, being the Minkowskian metric tensor for a D- 
dimensional spacetime. In (j2.ip . a is the AdS curvature radius which is related to the Ricci scalar by 
R = —D{D + l)/a^. In addition to the radial coordinate y we will use the conformal coordinate z, 
defined as z = ae^l"" . With this coordinate, the AdS line element is written in conformally-flat form: 

ds^ = {a/zf (r]ikdx'dx^ - dz^'^ . (2.2) 

The hypersurfaces z = and z = oo correspond to the AdS boundary and to the horizon, respectively. 
The dynamics of a fermionic field in curved spacetime are governed by the covariant Dirac equation 

il^^i,^ _ rnV' = , = 5^ + r^, (2.3) 

where F^ is the spin connection. The Dirac matrices 7^ are expressed in terms of the flat space- 
time gamma matrices 7^"\ as 7^ = e^^j7^"\ with e^*^^ being the tetrad fields obeying the relation 

^(a)^(6)^"'' ~ ' '^^'^ connection, one has 

= {l^^'h^'^el^^e^^,,, , (2.4) 

where the semicolon means covariant derivative of vector fields. For the geometry under consideration 
the tetrads can be taken in the form e^^^ = 6^z/a. With this choice, the spin connection has the 
following components (no summation over /) 

rD = 0, F; = g7(^)7«, l = 0,...,D-l. (2.5) 

For the combination appearing in the Dirac equation, we have 7'^r^ = —Dj^^^/{2a). 

We are interested in the change of the properties of the fermionic vacuum induced by the presence 
of the two boundaries, which are located at z = zi and z = Z2, z\ < Z2. The corresponding values of 
the physical radial coordinate, y, will be denoted by yi and y2- Zj = ae*/", j = 1, 2. We will assume 
that, on the boundaries, the field obeys the bag boundary condition, namely 

{l + ij^n'^^^)iP = 0, z = zj, (2.6) 

j = 1, 2, with n^^ being normal to the boundaries, = {—lyS^a/z. From these conditions it follows 
that the normal component of the fermion current vanishes at the boundaries. 
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In (D + l)-dimensional flat spacetime the Dirac matrices are N x N matrices, with N = 2^^^^^^^'^^ 
(the square brackets mean the integer part of the enclosed expression). In the discussion below we 
will assume the following representation for these matrices: 



with a = 1, . . . ,D and (Tadb + abd'a = '^^ab- The last relation directly follows from the anticommutation 
relations for the Dirac matrices. The matrices (j2.7p are related to the gamma matrices in the standard 



As we will see below, with the representation (j2.7p the equations for the components of the fermionic 
field are conveniently separated. This fact is to be remarked, since it allows for a complete calculation. 

Among the most important characteristics of the fermionic vacuum are the FC and the VEV of 
the energy-momentum tensor. For the evaluation of these quantities we will use the direct mode 
summation approach. In this approach we need a complete set of solutions to Eq. ()2.3p obeying the 
boundary conditions (j2.6p . For the positive-energy solutions the dependence of the mode functions 
on the time and on the coordinates parallel to the boundaries, denoted as x = (x^, . . . , x^~^), can be 
expressed in the form g*''^"*'^*, where k = (/ci, . . . ,/cd_i), kx = kix\ and the summation runs over 




(2.7) 



Dirac representation, by 7^"^ = iB^^y where 







(2.9) 



the initial Dirac equation is reduced to the set 




(2.10) 



from where we can obtain separate equations for the upper and lower components: 



[z'^dl - Dzd^ + X^z^ + D'^/A + D/2- m^a^ ± ar,ma) = 0, 




being 



up' — . With the substitution 




(2.12) 



Eq. (I2.11j) is reduced to the Bessel equation: 



( 



z^dl + zd^ + y?z^ ± aoma - rr?c? - 1/4) x± = 0. 




Taking the gd matrix in the form 





we further decompose X±(^) into upper and lower components. 




(2.15) 
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The solutions for these components dhectly fohow from (|2.11|) . and are 



(2.16) 



where Ju{x) and Yy{x) are Bessel and Neumann functions, respectively. The relation between the 
coefficients in these linear combinations are obtained by using Eq. (|2.10p . In order to find them we 
note that, from the anticommutation relations for the matrices ao and I = 1,...,-D — 1, and from 
(j2.14p . it readily follows that the matrices cj; have the form 



with bici = Cik = 1 and biCk = -bkCi, Cibk 
()2.10p with the upper sign, we find 



bi 
ci 



-Cf^bi, for / 7^ k. Using these relations, from the equation 



kkd^^^ + iXC_ 



t ' 



As a result, the solution of the Dirac equation can be written in the form 



AD+l)/2 



( 



\ 



UJ 



with the notations 
and 



\ hciZ^,u-i{Xz) +iXZ^^y_i{Xz) J 
V = ma + 1/2, 



c\.f Jt^iXz) + C^^Y^{Xz), 
c['lj^iXz) + d^^\iXz). 



Zi^Xz) 

We now must impose the boundary conditions (|2.6p . From the BC at z = zi, we have 



C 



+t 



c 



'+i _ Ju{Xzi) 



C 



(J) 



c 



(J) 



and hence 



Here and in what follows we use the notation 



Y,{Xzi)' 
= Z^{Xzi,Xz). 



Y,{x) 



(2.17) 



(2.18) 



(2.19) 



(2.20) 



(2.21) 



(2.22) 



In order to write the solution of the Dirac equation in a more compact form we introduce the 
notations 



{J) 



C 

c 



(J) 
+t 
{J) 

+i 



(2.23) 
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and 

With these notations, the solution ()2.17p obeying the boundary condition at z = zi can be expressed 
in the form 

- ^ ^(D+i)/2gikx~j<^t / Z_{Xzi,Xz)C^ \ ^2 25) 

\io-^Z+{Xzi,\z){iX + kiai)ci^^ J ' 

This solution corresponds to a state of the fermionic field with a given value of the momentum parallel 
to the boundary and with a given value of A. In order to completely specify the solutions we still 
need an additional quantum number. This corresponds to fixing the spinor Here we take 

= C^^^w^'^\ where C^^'' is a constant and w^'^\ a = 1, . . . , N/2, are one-column matrices of 

N/2 rows, with elements w^^^^ = Si^j. It can be seen that with this choice, the solutions (|2.25p . in 
combination with the corresponding negative-energy solutions (see below), do form a complete set 
specified by the quantum numbers (3 = (k, A, o"). Hence, the positive-energy mode functions obeying 
the boundary condition at z = zi have the form 



= 4+).^e^^™* ^ Z^{Xz, Xz)w^'^) y 

^ ^ ^ ^Z+{Xzi,Xz){iX + kiai)w^-) ' ^ ' 



Now we impose the boundary condition on the right boundary, located at z = Z2- Prom this 
condition it follows that the eigenvalues of the quantum number A are roots of the equation 

(7^,^_i(Azi,Az2) = 0, (2.27) 

where we have defined the function 

5.,^(x, y) = Ju{x)Y^{y) - J^,{y)Y,{x). (2.28) 

Eq. (I2.27P has an infinite number of positive roots. We will denote them by A = A„ = ^v,n/z\, 
n = l,2,.... 

The coefficient cl^'* in (I2.26P is determined from the normalization condition 



<p-^x / dz = <5(k - k')<5..'<5w, (2.29) 



22 

^ - / I --1/ I /I / r 

where 7 is the determinant of the spatial metric, I7I = {a/z)^ . By using a standard result for the 
integral involving the square of the cylinder functions (see, for instance, |d5j), for the normalization 
coefficient we find 

where we have introduced the notations 

ri = Z2/zi, (2.31) 

and 



Tu{r],x) 



This finishes the construction of the positive-energy mode functions. 

The negative-energy mode functions can be obtained in a similar way. They have the form: 



(2.32) 



^(-) ^ c'^-),^^^^-^^t ( ^^-(A^i, A.) {^X - ha,) \ ^^.33) 
P P \ Z^{Xzi,Xz)w^'^' j 

where is given by the same expression (j2.30p . As in the case of the positive-energy modes, for 

the eigenvalues of A one has A„ = '^v,n/z\- 
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3 Fermionic condensate 



In this section we consider the FC defined as the VEV (0| iptp |0) = {il^tp), where |0) corresponds to 
the vacuum state and ip = ip +7(0) is the Dirac adjoint. Note that the Dirac adjoint is defined through 
the flat spacetime matrix 7'^'^^ In addition to describing the physical structure of the quantum field 
at a given point, the FC plays an important role in models of dynamical chiral symmetry breaking 
(see the reviews [16] , for chiral symmetry breaking on curved spacetime with nontrivial topology, and 
|17j for recent developments). 

Expanding the field operator in terms of the complete set of positive- and negative-energy mode 
functions {ip^^\'>p^^ ^}, and using the anticommutation relations for the annihilation and creation 
operators, we find the mode-sum formula for the fermionic condensate: 

(VS^) = ^V;(-)4"\ (3.1) 

where 

„ N/2 00 

T.= dl.Y.E- (3-2) 

/3 a=ln=l 

The expression on the right-hand side (rhs) of (13. ip is divergent and, to make sense of it, some 
regularization procedure is needed. Here we assume that a cutoff function is present, without explicitly 
writing it. The special form of this function will not be important in the further discussion. 

Substituting the mode functions (12.330 into (13. ip after some transformations the fermionic con- 
densate can be expressed as 

(^/;^/;) = - 3 ^ / d\^^= 9v,vV1v,n,lv,nZ I z{)gy^y_x{lv,n,lv,nZ I Zx), (3.3) 

iOi^TTj Z^J n=l ^lln + k^zl 

the function gu,fiix, y) being defined by (|2.28p . As the roots -^y^n are given implicitly, the form (j3.3p for 
the FC is not convenient for the investigation of the effects induced by the boundaries. In addition, 
the terms in the series are highly oscillatory for large values of n. 

A more convenient form of the mode sum for the FC is obtained by using the summation formula 



{x,r]x) [h{ix) + h{—ix)] , (3-4) 



where 



"W(^.!/)= ,., ./^Y'"' (3.5) 

Ky{x)Iu-i{y) + Iu(x)Ky-i{y) 

and Iu{x), Ky{x) are modified Bessel functions. This formula is derived in Ref. [18j by using the 
generalized Abel-Plana formula (see also [E]). The corresponding conditions on the function h{u), 
analytic on the right half plane of the complex variable u, can be found in [18]. Applying expression 
()3.4p to the series over n in (j3.3p . after the integration over the angular part of k, and introducing a 
new integration variable x = kzi, the FC can be written as 

,r,> ,T,> SNiz/zif+'a-" , 0-2 

= <'»'^''+ (4,)(°^iW((p-l)/2) / 







00 



X / du — -===-Gy^y[u,uzl zx)Gy^y-x(u,uzl zx). (3.6) 

J X \ U — X 

where we have introduced the notation 

G,,^(x,y) = h{x)K^{y) - (-l)'^-%(y)if,(x). (3.7) 
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The first term on the rhs of (|3.6|) comes from the first integral in (|3.4p and it is given by 



{^ih)i = —- — — — — dxx du 

^ ' (47r)(^-i)/2r((D - l)/2) io Jo 

V? gu,u {u, uz/ zi )gu,u-i{u, u z/ zi ) ^ 



We now first consider this term. 



3.1 Condensate for the geometry of a single boundary 



The second term on the rhs of (13. 6p is finite for z < Z2, in the absence of the cutoff function. This 
term vanishes in the limit Z2 ^ oo, whereas does not depend on Z2. This allows us to interpret 

the part (13. Sp as the FC in the region z > zi for the geometry of a single boundary at z = zi when the 
other boundary is absent. This can also be seen by direct evaluation of the FC using Eq. (|3.ip . The 
corresponding positive-energy mode functions are given by (j2.26p . where now the spectrum for A is 
continuous. Consequently, in (13. ip we have = I lo^ ^^Yla=i- '^^^ normalization coefficient is 
determined by the condition which is obtained from (j2.29p by the replacements f^^ dz — )• j"^ dz and 
^nn' — 5- (5(A — A'). In this way, we can see that 



^ ' 2(27r)^-i 



1 + 



Y^{\z^] 



(3.9) 



Similarly, the negative-energy modes have the form (j2.33p with given by the same expression 

()3.9p . Substituting these mode functions into the formula (j3.ip . the expression (|3.8p is obtained for 
the single boundary part. 

For further transformation of such expression we use the identity 

9uA'^^y)9u,u-i{u,y) T f \ T / \ ^ Sr- •^'^(^) u(s)^ MjW / \ ^Q1n^ 

+ = ■'-fa>''-'(!'> - 5 E Hy)E^-M, (3.10) 

the Hjf\y), s = 1,2, being Hankel functions. Substituting (j3.10p into (j3.8p . in the integral over u 
with the second term on the rhs of (j3.10p we rotate the integration contour in the complex plane u by 
the angle tt/2 (— 7r/2) for the term with s = 1 (s = 2). Introducing the modified Bessel functions, the 
FC is then decomposed into 

{i;iJ), = {^PiJ)o + {iJi^)f'\ (3.11) 

with 
and 

/7/\(b) N{z/zr)^+-'a-D . D-2 

X- / du^==-p-^^K^{uz/zi)K^-i{uz/zi). (3.13) 

The second part (j3.13p is finite for z > zi and the cutoff function need be kept in the first part (|3.12p 
only. Note that the latter does not depend on z. The term {'i/jip)^^^ vanishes in the limit z — )■ oo 
and, hence, the term {'4>ip)o can be interpreted as the pure AdS part of the FC when the boundaries 
are absent. The property that this term is uniform is just a consequence of the maximal symmetry 
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both of AdS spacetime and of the vacuum state we have chosen. With the representation (|3.1ip . the 
renormahzation of the FC outside the boundary is thus reduced to the one corresponding to AdS 
spacetime when the boundaries are absent. 

The expression (j3.13p for the boundary-induced part can be further simphfied by using the inte- 
gration formula 





dxx^-^ di, ' = — / drr"-^ fir). (3.14) 

' V^^^ ^ 2Y{D/2) Jo ■'^ ' ^ ' 



This is obtained by introducing a new integration variable y = \Jv? — x^, passing to polar coordinates 
on the plane (rc, y) and then integrating over the polar angle. Now, the expression for the FC induced 
in the region z > Zj by a single boundary, located at z = zj, simply reads 

(b) ^ _2Nizlz^ r dxx^^KAxz/z,)K,^^{xz/zj), (3.15) 
where we have introduced the notation 

Ad = {Ai:)^''^T{D/2). (3.16) 

Note that the boundary-induced part (|3.15|) is a function of the ratio z/zj alone. By taking into 
account that z/ Zj = e^y~y3^l°- ^ we see that for a given distance from the boundary, y — yj, the quantity 

{'tpip)f'^ does not depend on the location of the boundary. Again, this is a consequence of the maximal 
symmetry of the AdS bulk. It follows from (j3.15p that (Tpi/j)^^'^ is negative. 

The expression (j3.15p gives the boundary- induced part in the region z > Zj. In order to find the 
condensate induced by a single boundary in the region to the left of the boundary, we take the limit 

— )• in the general expression (j3.6p for the geometry with two boundaries. Introducing in (j3.15p . 
with j = 1, a new integration variable, y = x/zi, we see that in this limit {ijjip)f'^ vanishes and (ipip)! 
reduces to (V'V')o- Iii the second term on the rhs of (13. 6p . passing to the new integration variables 
u' = u/zi and x' = x/zi, the limit zi — )■ is readily evaluated. As a result, from p.6p we obtain the 
FC in the region z < Z2 for the geometry of a single boundary at z = Z2- For the geometry of a single 
boundary at z = zj, the FC in the region z < zj is expressed as {'ipip)j = (V'V')o + {'^''P)f'\ with the 
boundary-induced part being 

2N{z/z,)^+^ r^^^^^ K^^.jx) 



= -tliy^pj^ I dxx^tlI^^U{xz/z,)U^,{xz/z,). (3.17) 



This expression is finite for points away from the boundary and vanishes at the AdS boundary. Similar 
to the case of the region to the right, the FC defined by (j3.17p does not depend on the location of the 
boundary for a fixed distance from it. Note that the FC is not symmetric with respect to the boundary, 
the reason for this being that, though the background spacetime is homogeneous, the boundary z = Zj 
has a nonzero extrinsic curvature tensor and the two sides of the boundary are not equivalent. Here 
the situation is similar to that for curved boundaries in the Minkowski bulk. 

Let us consider the asymptotic behavior of the single boundary-induced part. At large distances 
from the boundary, z ^ Zj, we introduce in (j3.15p a new integration variable y = xz/zi and expand 
the integrand by using the formulas for the modified Bessel functions for small values of the argument. 
To leading order, the remaining integral which involves the product of two MacDonald functions is 
evaluated by using a formula from [15) . In this way, one finds that 

/,Wa(^) ~ Na-^{z,/zf^ DT{D/2 + 2u)nD/2 + v) 

\'P'P)j - 2^+2^+i7r(^-i)/2 uT{D/2 + iy + l/2)T^u) ' ^ ' 



9 



For the region to the left of the boundary, assuming z <C by direct expansion of the integrand in 
()3.17p . to leading order one gets 



This expression gives the asymptotic behavior of the FC near the AdS boundary. Note that both 
limits, ()3.18p and ()3.19p . correspond to distances from the boundary much larger than the curvature 
scale for the background spacetime: \y — yj\ ^ a. As we see, in these regions the boundary-induced 
part decays exponentially with the distance from the boundary. 

For points near the boundary, \\ — zj Zj\ <C 1, the dominant contributions to the integrals in (j3.15p 
and ()3.17p come from large values of x. By using the asymptotic expressions for the modified Bessel 
functions for large values of the argument, to leading order we find 

This leading term coincides with the corresponding one for the boundary in Minkowski spacetime. 
3.2 FC in the region between the two boundaries 

Now we return to the geometry with two boundaries. Using of the integration formulas (|3.14p in (j3.6p . 
the FC in the region z\ < z < Z2 can be expressed in the form 



Aoa^ JO 

xQ.^J-\x,r]x)G,y^i,{x,xz/zi)G^^,y^i{x,xz/zi), (3.21) 

where {'ipip)^^'^ is given by (|3.15p with j = 1. For points outside the boundaries renormalization is 
only needed for the first term on the rhs. In formula ()3.2ip . the functions Q,ly\x,y) and Gy^y^i{x,y) 
are positive. Using the result that the function Iy{x) / Ky{x) is a monotonically increasing in terms 
of X > 0, we see that the function Gu,u{x,y) is negative for x < y and positive for x > y. Hence, 
the last term in (j3.2ip is negative. Combining this with the result {xp'ip)^''^ < 0, we conclude that the 
boundary-induced part of the FC in the region between the boundaries is negative. 

It can be checked that, in the region between the boundaries, the FC can also be written as 



with 







xn^^\x/r], x)Gu-i,u-iix, xz/z2)Gu-i,y{x, xzjz^), (3.22) 



'^^'V)' ^ , J'lfr^fl^ ,y (3.23) 

Ky[x)Iv-i{y) + Iu[x)Ky-i{y) 

In (|3.22p . {tjjTp)^^ is given by p.l7p with j = 2 and the last term is induced by the boundary at 
z = zi. The latter vanishes in the limit zi — ?• 0, it is finite for z = Z2, and diverges on the left 
boundary. Divergences occurring in this term are the same as those for a single boundary at z = Z2- 
All functions in the integrand of the last term in (|3.22p are positive. From the discussion above it 
follows that, for the region between the boundaries, if we express the FC in the form 

.7 = 1,2 
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then the last (interference) term is finite everywhere, including the points on the boundaries. The sur- 
face divergences are contained in the single boundary parts. At large distances between the boundaries, 
as compared with the curvature radius of the background spacetime, (?/2 — Vi) ^ a, the interference 
part is exponentially suppressed. 

Consider now the Minkowskian limit of the expression for the FC in the region between the 
boundaries. This corresponds to ma ^ 1 for a fixed value of y, hence, one has z/ ^ 1 and z/zj ~ 
1 + [y — yj) / a. Introducing in the formulae above a new integration variable u = x/v and using the 
uniform asymptotic expansions for the modified Bessel functions for large values of the order, after 
some transformations, to leading order we get 

±j J m x—m 

x{(m + x)[e2^(2^-S'i) +e2^(^2-y)] _ 2m], (3.25) 

where (V'V')(m) is the FC for boundaries in a Minkowski bulk. The expression (|3.'25p for (V'V')(m) is 
a special case of a more general formula derived in [20] for the Minkowski bulk with compact spatial 
dimensions. The fermionic condensate for a massless field has been considered in j21j . 



4 VEV of the energy- momentum tensor 

The VEV of the energy-momentum tensor is another important local characteristic of the fermionic 
vacuum. In order to find this VEV we use the mode-sum formula 

. „ N/2 oo 
(T=l n=l 

where the brackets denote symmetrization over the indices enclosed. Note that for the covariant 
derivative of the Dirac adjoint field one has V^V^^ ^ = i9^V'^ ^ — V'^ ^F^. We see that the spin connection 
appears in the expression for the VEV in the form ^[^^u) + '^{^ilv)- By using the expressions ()2.5p it 
can be seen that this combination vanishes. Hence, in the evaluation of the VEVs we can make the 
replacement V y — )• dy. 

Substituting (I2.33P for the negative-energy mode functions into (14. ip . it can be seen that the 
off-diagonal components vanish. For the VEVs of the diagonal components we find the following 
expressions (no summation over /i) 

^^^^"4(47r)(^-i)/2r((Z)-l)/2)yo ^1 .fii^^ 

with the notations 

f!)°Hx,y) = - {x^ + u^) [9l,uix,y) + gl^^_i{x,y)] , 

f^'\x,y) = -^[4^(a;,y)+4,_i(x,y)], (4.3) 

/^^Hx,y) = 



9lJx,y) + 9l,^-i{x,y) - — — -guAx^y)9u,u-i{x,y) 

. y 



and / = 1, . . . , D — 1. As in the case of the FC, here we assume the presence of a cutoff function which 
makes the expressions (|4.2p finite. 

Further transformations of the VEVs proceeds similarly to those in the case of the FC. First, the 
series over n in (j4.2p is transformed via the summation formula (j3.4p . Next, for the part corresponding 
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to the second term on the rhs of (|3.4p we use the integration formula (|3.14p . In this way, the VEVs 
are expressed as (no summation over fi) 



N{z/zi) 

where the notations 



D+2 foo 



TO = TOi + ^5^ J dxx^+'ni'\x,rjx)Fi^\x,xz/z,), (4.4) 



FS{x,y) = :i[Gl^ix,y)-Gl,_,{x,y)], 

F[^\x,y) = G2^_i(x,y)-G2^(x,y) + ^^G,,,(x,y)G,,,_i(:r,y), (4.5) 



are introduced with / = 0, . . . ,D — 1. In (14.41) . the first term comes from the first integral in the 
summation formula ()3.4p and it has the form (no summation over ^) 



^ /i^^(x,xz/zi) 
^^V^2^j2(x) + y,2(^)- ^^-^^ 

This term corresponds to the VEV of the energy-momentum tensor in the region z > z\ for the 
geometry of a single boundary at z = zi. The second term on the rhs of (14.41) is induced by the 
presence of the second boundary at z = Z2. The latter is finite for z\< z < zi and renormalization is 
needed for the first term, only. 

The single boundary part (|4.6p could also be directly obtained by using the the corresponding 
mode functions. The latter are given by ()2.26p and (I2.33p . with < A < and with the normalization 
coefficient (13. 9p . 

The transformation of the part (T^)i is also similar to that for corresponding term in the FC. 
Using ()3.10p and the similar identity 

9l,^{x,y) 1 ^ J^(x) 



J; 



Jl{x) + Y^{x) ^tt2^^^^(x) 

the integrand in (j4.6p can be decomposed into parts containing Bessel and Hankel functions, with the 
argument xz/zi. Then, rotating the integration contour over x by an angle 7r/2 (— 7r/2), for the part 
with the function H^\xz/ zi) {Hj^\xz/zi)), we get (no summation over //) 

{TPi = {TPo + {Tpf'\ (4.8) 
where the two terms are respectively given by 

2r((Z) - l)/2)a^+i io Jo V^^T^' ^ ' 

and 



W>1 --^l^X </..-'^S!-(../.,). (4.10) 

In ()4.9p . the expressions for fQ^{y) are obtained from the corresponding expressions for fl)^\x,y) by 
the replacement gu,a{x,y) — )• Ja{y), and in ()4.10p we have defined 

S{2i^) = ^ [K'^ix) - K^,_,{x)] , 

S[^\x) = Kl_^{x)-Kl{x) + '^^^K,{x)K,.i{x), (4.11) 
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with / = 0, . . . , D — 1. As we see, in both cases of single and double boundary geometries the energy 
density is equal to the stresses along the directions parallel to the boundaries. This property is related 
to the symmetry of the problem under consideration. 

ya^j^ishes in the limit zi — )• (for the corresponding asymptotic behavior see below) while 
(Tq )o can be identified with the VEV of the energy-momentum tensor in the boundary-free AdS bulk. 
With the decomposition (j4.8p and for points away from the boundary, renormalization is required for 
this last part only. Because of the maximal symmetry of the background geometry, the renormalized 
VEV does not depend on the spacetime position and is completely determined by the trace: (Ta)o = 
{T^)o6a/{D + 1). For the case D = 3 this VEV was investigated in [22j using the zeta-function 
techniques and also Pauli-Villars regularization. In what follows we specifically discuss the effects 
induced by the boundaries. 

For a massless field one has u = 1/2 and from (j4.10p it can be easily seen that (r^)^''^ = for 
z > zi. We could obtain this result directly, by taking into account that for a massless fermionic 
field the problem is conformally related to the corresponding problem for a single boundary in the 
Minkowski bulk and for the latter geometry the VEV of the energy-momentum tensor vanishes. In the 
region between two boundaries, by making use of the expressions G1/2, 1/2(2^) v) = {xy)~^^'^ sinh(j; — y) 
and ^1/2,-1/2(2^)^) = (xy)~^/^ cosh(a; — y), from (14.41) we get (no summation over /x) 

{TP = {Tli)o + {^y^' (^/f)(M),„^=0. (4.12) 

where 

{TP(M),m=o - - (4^)(D+i)/2(^2- ^1)^+1 ^""^ ^ ' 

for fi = 0, . . . , D — 1, and {Tl^)(^M),m=o = "-^(^o )(M),m=o- (|4.12p . Cr{^) is the Riemann zeta 
function. Here, {Ta)(M),m=o is the corresponding VEV for two boundaries at z = zi and z = Z2 in 
Minkowski spacetime. Of course, (j4.12p shows the standard conformal relation between the problems 
in AdS and Minkowski bulks. Note that for a massless field the boundary-free part (Tq )o is completely 
determined by the trace anomaly. 

We can write the VEV of the energy-momentum tensor in an alternative form (no summation over 

{TP = {TPo + (T^n?^ + ^Aoa^T dx x''^^n^\xh,x)Fi^ {x,xz/ z^), (4.14) 
where the term 

is the part in the VEV induced in the region z < ^2 by a single boundary at z = Z2 when the boundary 
at 2; = zi is absent. In (j4.14p and (j4.15p we have introduced the notations 

Si^\x) = l2(x)-/2_,(x) + ^^/,(x)/,-i(x), (4.16) 



and 



T^i^^y) = [Gl-i,u^iix,y) -Gl_^„{x,y)] , 
T2u\x,y) = Gl^i„{x,y) - Gl_i„_i{x,y) - '^^ ^ ^ Gu-iM^i{x,y)Gu-^i^^{x,y), (4.17) 

with I = 0, . . . , D — 1. The last term in ()4.14p is induced by the right boundary. It is finite for 
zi < z < Z2 and diverges at z = zi. It can be seen that, for a massless field one gets (no summation 
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over fx) {Tli)^2^ = (2;/a)^^^(T^)(M),m=0! where {Tli,){M),m=o is given by (|4.12p with zi = 0. As we see, 
in the region z < Z2 the problem with a single boundary in the AdS bulk is conformally related to 
the problem in the Minkowski spacetime with two boundaries. This is a consequence of the boundary 
condition we have imposed on the AdS boundary. 

It can be checked that both the single-boundary induced part, {Ta)^^\ and the second boundary- 
induced part (last terms in (I4.14j) and M.lSh ) in the VEV of the energy-momentum tensor obey the 
trace relation = mipip. In addition, the VEV obeys the covariant continuity equation Ta-^ = 
which, for the geometry under consideration, reduces to the single equation 

z^+^d, {tE/z^+^) +TI^ = 0. (4.18) 

Using the trace relation and taking into account that the boundary-induced part in the FC is negative 
everywhere, from (|4.18|) we conclude that the boundary-induced part in the VEV of the normal stress 
is a monotonically increasing function of z for all points outside the boundaries. 

Taking into account that Ii^{x) < Iy_i{x) and Ky{x) > K^_i{x), from ()4.10p and ()4.15p we see 
that, for the geometry of a single boundary at z = zj, the boundary-induced part in the energy density 
is negative everywhere: {Tq)^^^ < 0. Using the properties of the modified Bessel functions, it can be 
seen that Sl^\x) > for z > Zj and Sj^\x) < for z < Zj. From here it follows that (T^)^'^^ < in 
the region z > zj and (Tj^)^^^ > in the region z < zj. Next, it can be checked that 

Gu,,y^i{x,y) > -Gy^y{x,y) > X <y, 

Gu-.i^y{x,y) > Gu-i,u-i{x,y) > 0, X > y, (4.19) 

From these relations we see that Fj'^\x,y) < and the parts of the energy density induced by the 
second plate (last terms in (|4.4p and (j4.14p with fi = 0) are negative. Hence, for the geometry of two 
boundaries the energy density is negative everywhere. 

Now let us consider the asymptotics for the single boundary parts in the VEV of the energy- 
momentum tensor at small and large distances from the boundary. For the boundary at z = Zj, at 
large distances, z ^ Zj, to leading order one has 

,^Dab) ^ DiTS)f^ ^ Nma~^iz,/zr r{D/2 + z.)r(g/2 + 2^.) 

^ D + 2u 2^+2'^+%(^-i)/2i/r2(zy)r(L»/2 + z/ + 3/2)' ^ ' 

with Zjl z = e~^'^~y^^^°' . For the region z < Zj with the condition z Zj, the leading order terms have 
the form: 

^^^^^ ^-2^^^o^- ^ 22-M^.r2(.) X TZM- ^ ^ ^ 

The relations between the energy density and the normal stress, given by (j4.20p and (j4.2ip . can also 
be obtained by using the continuity equation (lilHl) for (r^)^''^ with {Tll:)f^ = D{T^)f^ + {Tg)f\ 
As it is seen from (j4.20p and (j4.2ip . at distances from the boundary larger than the AdS curvature 
scale, the boundary-induced part in the VEV of the energy-momentum tensor decays exponentially. 
For a scalar field with curvature coupling parameter ^ and with Robin boundary condition at z = zj, 

for the VEV of the energy-momentum tensor one has [TO]: (7/^)^-^'' oc {zj / z)'^'^^" for z ^ Zj and {Tji)^^^ 
oc (2;j/2;)^+2j^sc ^qj. 2 ^ 2^.^ where u^c = \/ — D{D + 1)^ + rn^a?. In particular, for a conformally 
coupled scalar one has Usc = ■\/m?a?' + 1/2, and the suppression of the VEVs is weaker than in the 
fermionic case with the same value of the nonzero mass. 

For points near the boundary, the dominant contributions to the integrals in (|4.10p and (j4.15p 
come from large values of x. By using the asymptotic expressions for the modified Bessel functions 
for large values of the argument, we get the following leading behavior: 

0^ W ^ {D-^)a ,rj.D,{b) ^ Nmn{D + l)/2) 
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For a boundary in the Minkowski spacetime the leading terms for the energy density and stresses 
along directions parallel to the boundary coincide with (j4.22p . whereas the normal stress vanishes. 

The Minkowskian limit of the formulas for the VEV of the energy-momentum tensor is taken in 
a way similar to that we used for the case of FC. For ma ^> 1 and for a fixed value of y, to leading 
order we find (no summation over //): 

{TP - (T,-)(M) = I dx ^^,J_,,,^^ G[%ix), (4.23) 

where 

for /X = 0, 1, . . . , L> - 1, and Gj2j(x) = -2x^. Here, (T^) (M) is the VEV for the geometry of two 
boundaries in the Minkowski bulk. The formula (|4.23|) is a special case of the result given in |2U| . 
The fermion Casimir energy for two parallel plates in 4-dimensional Minkowski spacetime has been 
investigated in [23j and [24] for massless and massive fields, respectively. The corresponding result for 
arbitrary number of dimensions is generalized in ^25j. The topological Casimir effect and the VEV 
of the fermionic current for a massive fermionic field in a spacetime with an arbitrary number of 
toroidally compact dimensions have been considered in |26j . 

In Fig. [H for the geometry of a single boundary, located at y = 0, we display the boundary-induced 
parts of the VEV of the energy density (n = 0, full curve) and of the normal stress (// = D, dashed 
curve), as functions of y/a. The latter measures the distance from the boundary in units of the AdS 
curvature radius. The graphs are plotted for a fermionic field in 4-dimensional AdS spacetime {D = 3) 
and for the mass we have taken ma = 1. 
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Figure 1: Boundary- induced parts of the VEV of the energy density (full curve) and of the normal 
stress (dashed curve), induced by a single boundary, at y = 0, as functions of y/a. Both plots are for 
D = 3 and ma = 1. 



5 Interaction forces and the Casimir energy 

The force acting per unit surface of the boundary at z = Zj (vacuum effective pressure) can be obtained 
evaluating the normal stress at the location of the boundary: p^^') = — l^T£)z=Zj- The boundary-free 
parts of the force acting from the left- and from the right-hand sides of the boundary compensate and 



y) 



(4.24) 
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the resulting force is determined by the boundary induced part. By using the decomposition for the 
VEV of the energy-momentum tensor, in the region between the boundaries, the effective pressure is 
obtained as 

(5.1) 

where the first term on the rhs correspond to the situation when the second boundary is absent and 
the second term is induced by the presence of the second boundary. The latter can be termed as the 
interaction contribution. For the first part, one has pf^ = —{T£)j^z=Zj- In the regions z < zi and 

z > Z2, the single boundary terms remain only: p^^^ = p^^^ for z < zi and p^'^^ = p^^^ for z > Z2- 

(i) 

Because of the surface divergences in the single-boundary parts of the VEV, the term p^ is divergent 
and needs renormalization. The interaction part is finite for all nonzero values of the distance between 
the boundaries and it is not affected by the renormalization procedure. 

The interaction parts of the effective pressure are obtained from the last terms in (j4.4p and (I4.14p . 
by setting in the expressions for the |^-components z = zi and z = Z2, respectively. Using the 
Wronskian relation for the modified Bessel functions, one finds 

As we see, the corresponding effective pressures are always negative and, hence, the interaction forces 
are attractive. Using now the relations 



^ ^ Iu{xzi)Ky^i{xz2) 

Ky{xZl)ly^l{xZ2) 

the interaction parts can be expressed in a more convenient way as 



(i) _ N{zj/a) 



D+l 



P{int) 



{-IVAd 



I dxx^ ^ In 





1 + 



Iu{xzi)K„^i{xz2) 

Ku{xZi)Iu-l{xZ2) 



(5.3) 



(5.4) 



Note that the forces acting on the left and on the right boundaries are different, in general. This 
property is related to the nonzero extrinsic curvature tensor for the the boundary geometry under 
consideration. 

At small distances between the boundaries, r/ — 1 ^ 1, the dominant contribution comes from large 
values of x and, to leading order, we get 

(,) NDT{{D + l)/2) 

^(int) - (47r)(^+l)/2(y2 - yi)^+^ ' ^ ' 

At large distances one has ^> 1, and the leading terms in the asymptotic expansions are given by 
the expressions 

(1) ^ D+2u^iKu-i{x)_ 

P(int) - ^^r2(^)^D+2. /.-l(x)' 

Now, let us consider the Minkowskian limit corresponding to ma ^ 1 for a fixed value of y. In 
this limit one has ~ 1 + (1/2 — Vij/o,- Introducing in (|5.2|) a new integration variable, u = x/v, we 
use now the uniform asymptotic expansions for the modified Bessel functions for large values of the 
order. After some transformations, to leading order we have 

(i) D+l r , X2(x2- 1)^/2-1 
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Of course, in the Minkowskian limit the forces are the same for the boundaries at y = yi and y = y2- 
Note that, for the geometry of a single boundary, in the Minowskian limit the stresses on the left- and 
right-hand sides are the same by the symmetry of the problem. As a result, the corresponding net 
force vanishes and only the interaction part remains. 

In Fig. [21 we depict the interaction forces for the model with D = 3 as functions of the separation 
between boundaries as measured in units of the Compton wavelength of the fermionic particle. The 
dashed curve is the corresponding force for boundaries in Minkowski spacetime. The plots are for 
ma = 0.5 (black), ma = 1 (blue) and ma = 2 (red). The curves on the left (right) of the dashed curve 
are for j = 1 (j = 2). As we see, with increasing a the forces tend to the corresponding result for the 
Minkowski bulk. 




0.0 0.5 1.0 1.5 

m{y2-yi) 



Figure 2: Interaction forces per unit surface as functions of the separation between the boundaries 
for the fermionic field in a 4-dimensional AdS bulk. The dashed curve corresponds to the force for 
boundaries in the Minkowski bulk. The plots correspond to ma = 0.5 (black), ma = 1 (blue) and 
ma = 2 (red), respectively. 

Now we consider the total vacuum energy in the region between the boundaries. The formal 
expression for this energy is obtained by integration of the energy density given by ()4.2p . with /i = 0. 
Making use of standard formula for the integrals involving the square of a cylinder function, we get 

'21 



N f dk 



2 J (2^)^-1 



E \^<n/zl + k^. (5.^ 



=1 



This formula expresses the vacuum energy as a sum of ground state energies for elementary oscillators. 
Obviously, expression ()5.8p is divergent and regularization, with the subsequent renormalization, is 
necessary. Here we follow the zeta function approach (for the application of the zeta function tech- 
niques to the calculations of the Casimir energy see [27] and references therein). Instead of (j5.8p . we 
consider the related zeta function 

Cis) = -Mo^4 / E {llJz! + kr'^' , (5.9) 

where the constant fiQ with dimensions of mass is introduced by dimensional reasons (it is the regu- 
larization parameter). The expression on the rhs of (|5.9p is finite for Res > D and for the evaluation 
of the vacuum energy we need its analytic continuation at s = —1: Ezi<z<z2 = C(s)|s=-i- 
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Evaluating the integral over k, the zeta function can be expressed in the form 



NTjis + l- D)/2)fi; 
'2(4^)(^-i)/2r(s/2). 



s+l 



n- 



— Ci (s + l-D), 



(5.10) 



where we have introduced the partial zeta function 

oo 



(5.11) 



n=l 



directly related to the eigenvalues 71,^,1. We need the analytic continuation of the function ()5.1ip to a 
neighborhood of s = —D. This procedure is standard in the theory of the Casimir effect and we only 
give the main steps. 

From Cauchy's residue formula, the integral representation follows 



Jc 2v^^ 



(5.12) 



where C is a closed, counterclockwise contour on the complex z plane enclosing all zeros 'yu,n- We 
take the contour made of a large semicircle (with radius tending to infinity) centered at the origin and 
placed to its right, plus a straight part overlapping the imaginary axis and avoiding the origin by a 
small semicircle Cp on the right half-plane with radius p. For small p, one has 



^ u In [ug^^_i{u, ur])] = 



■ sin(7rs/2), 



(5.13) 



being B a constant independent of s. Denoting the upper and lower halves of the contour C by Ci 
and C2, respectively, the integral can be cast in the form 



Cl is) 



du 



Ir, 27ri 



=1,2 



Cl 



2TTi 



u ^9uln 



j.(^)ga(^^) 

Hu°'\u)Jy-i{u'qz) 



(5.14) 



After parameterizing the integrals over the imaginary axis and substituting into (j5.10p . we arrive at 
the following expression for the zeta function 



C(^) 



(4vr: 



l(l-D)/2 



+l^l+s--D 



2T{s/2)T{{D + 1- s)/2) \D + l-s 



irBp 



D+l-s 



+ 



dx X 



+ ln(x'^i^^(x))+ln 1 + 



In {x^-'' h^i{xj])) 

Iv{x) Kv-i{xvi) \ 



K{x)Iy-i{xr]) J 



(5.15) 



where we have used that T{y) sin(7ry) = 7r/r(l — y). The term with the factor B vanishes in the limit 
p — )• at the physical point s = —1, while the last term is finite at this point. The first (second) term 
in the square brackets in (I5.15P corresponds to the vacuum energy in the region z < Z2 {z > zi) for 
the geometry of a single boundary at z = Z2 {z = zi) and will be denoted as E^']^^^^ (-^iWzi)- Adding 
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also the vacuum energies from the regions z < zi and z > Z2 (denoted as e[^1<zi ^^'^ ^'\\>Z2)^ 
total vacuum energy we finally get 

(7) 

Here, -Bj is the energy for the geometry of a single boundary at z = Zj, and the interference part 
AE is given by the last term in (jS.lSp with s = —1. After integration by parts, one gets for the latter 

(5.17) 



^ _^ Iu{xzi)Ki,^i{xz2) 



K{xzi)Iy^i{xz2) 



The renormalization procedure to be carried out for the divergences of the single boundary parts in 
(|5.16|) is similar to those previously discussed within the framework of braneworld scenarios (see, e.g., 
O [8]). Note that the dependence of on Zj is in the form Zj ^ and this dictates the form of 
the counterterms located on the boundaries. By finite renormalizations, the single boundary terms 
in (|5.16p can be absorbed into the counterterms. The renormalized vacuum energy has the form 
E = Ylj=i 2 '^i^J^ ~^ AE, with renormalized coefficients cj. Comparing (I5.17P with (15. 4p . we find the 
relation 

between the interaction forces and the interference part of the Casimir energy. 

Here we have considered the fermionic Casimir effect with bag boundary conditions on the back- 
ground of AdS spacetime. In a similar way the fermionic Casimir densities can be evaluated for the 
Randall-Sundrum braneworld model. This model is formulated on 5-dimensional AdS spacetime, thus 
with a single extra dimension. The fifth dimension y is compactified on an orbifold /Z2 of length 
L, with —L < y < L. The corresponding line element has the form ()2.ip with the warp factor e"'^'^'/'^. 
Two 3-branes are located at the orbifold fixed points, y = and y = L. In terms of the conformal 
radial coordinate z, for the branes one has z = a and z = zl = ae^^"". In the Randall-Sundrum model, 
depending on the parity of the spinor field under a chiral transformation, two types of boundary con- 
ditions arise on the branes. For these boundary conditions the eigenvalues of the quantum number A 
are roots of the equation 

gu~s,u-s{>^a, Xzl) = 0, (5.19) 

with s = 0, 1 for even and odd fields, respectively. The summation formula for series over these eigen- 
values is obtained from the general results of |18l[l9]. Calculations are actually the same as those we 
have described above for the case of bag boundary conditions. In the normalization condition for the 
mode functions the integration over y goes over the region {—L,L). As a result the normalization co- 
efficient will have an additional factor 1/2, as compared with the case where the problem is formulated 
on the interval (0, L). 

The FC and the VEV of the energy-momentum tensor for a fermionic field in the Randall-Sundrum 
model are obtained from the formulas given in Sects. [3] and H] by changing the order of the appropriately 
modified Bessel function from v — 1 to v for even fields and from v to v — 1 for odd fields and by 
adding an extra factor of 1/2. For example, in the case of even fields, to Eq. (j3.15p we need only add 
the factor 1/2, whereas to Eq. (I3.17P we have to add the factor 1/2 and make also the replacements 
Kjj-i{x) — )• —Ku{x) and Iu-i{x) — Iu{x). In the case of odd fields the situation is just opposite: in 
(13.15P we add the factor 1/2 and replace K^{x) — )■ —Ky-i{x) and Iu{x) Iy-i{x), while in (13.17P we 
only add the factor 1/2. Note that when evaluating the vacuum energy in the braneworld model the 
integration goes over the region {—L, L) and there is no need to add the factor 1/2 in the corresponding 
expressions. 

Through the above mentioned replacements of the modified Bessel functions, from (j5.17p we readily 
obtain the corresponding vacuum energies for even and odd fields in the Randall-Sundrum model. For 
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D = 4 the corresponding formulas were obtained in |8J (note that in this reference the effective 
Lagrangian is considered per fermionic degree of freedom, which corresponds to —AE/N). The VEV 
of the energy-momentum tensor for a bulk Dirac spinor in the Randall-Sundrum model has been 
considered in [IT]. In this reference, for the case of a massive field, a general formula is given for 
the unrenormalized VEV only. To compare, in our approach, based on the generalized Abel-Plana 
formula, the pure AdS parts in the VEVs are extracted explicitly and, for the points away from the 
branes, the renormalization procedure is the same as for the boundary-free parts. In addition, the 
boundary induced parts are presented in terms of exponentially convergent integrals, which are very 
well suited for numerical calculations. 

6 Conclusions 

In this paper we have investigated the fermionic condensate and the VEV of the energy-momentum 
tensor for a massive fermionic field in AdS spacetime in the presence of two boundaries on which 
the field obeys bag boundary conditions. For the evaluation of the VEVs we have employed the 
mode summation technique. In the region between the boundaries, a complete set of positive- and 
negative-energy mode functions is given by (j2.26p and (j2.33p . respectively, where the eigenvalues of 
the radial quantum number A are determined from the boundary conditions and they are solutions of 
the equation (I2.27p . The mode sums for the FC and the energy-momentum tensor contain series over 
these eigenvalues. For the summation of the series we have used the generalized Abel-Plana formula 
()3.4p . which allowed to separate the VEVs into single boundary and second boundary-induced parts. 
In this representation, explicit knowledge of the eigenvalues of A is not necessary. The VEVs for 
the geometry of a single boundary are further decomposed into boundary-free and boundary-induced 
parts. As a result, in the region between the boundaries, the VEVs can be expressed in two equivalent 
ways, respectively given by Eqs. (I3.21j) and (I3.22p . for the FC, and by Eqs. (14. 4p and (I4.14p . for the 
energy-momentum tensor. With these representations, and for points away from the boundaries, the 
boundary induced part is finite and renormalization is required for the boundary-free part, only. 

For the geometry of a single boundary located at z = Zj, the boundary-induced contribution to 
the FC is given by (I3.15p . in the region z > Zj, and by (I3.17p . in the region z < Zj. This contribution 
is negative for both regions and it is not symmetric with respect to the boundary. Such fact is related 
to the existence of a nonzero extrinsic curvature tensor for the boundary z = Zj in the AdS bulk. 
At large distances from the boundary, as compared with the value of the AdS curvature radius, the 
boundary-induced parts are exponentially suppressed by the factors e^^^^^"^-'-'/'^, in the region y > yj, 
and e"^'^"'"^'^^^^"^-'-'/", in the region y < yj. In particular, the boundary-induced part vanishes on the 
AdS boundary. For points near the boundary at y = yj, the leading term in the asymptotic expansion 
of the FC is given by (j3.20p and it does coincide with the corresponding expression for the boundary 
in Minkowski spacetime. For the geometry of two boundaries, the FC in the region between them can 
be expressed in the form (j3.24p . where the interference term is finite everywhere, including the points 
on the boundaries. For large separation of the boundaries, as compared with the curvature radius of 
the background spacetime, the interference part is exponentially suppressed. The boundary-induced 
part of the FC in the region between the two boundaries is negative. 

The boundary-induced contributions in the VEV of the energy-momentum tensor for a single 
boundary are given by ()4.10p and (14.150 . for the regions on the right and on the left of the boundary, 
respectively. The corresponding vacuum energy is negative, whereas the normal stress is negative on 
the right domain and positive on the left one. At large distances from a single boundary located at 
y = yj, the boundary-induced terms decay as e"^'^^^"^-')/", for y > yj, and as e~^^~^'^^'>^y~y^'^/°'^ for 
y < yj. For points near the boundary the corresponding asymptotic behavior is given by (I4.22p . 

The forces acting on the boundaries and the Casimir energy were considered in Sect. [5l The 
force acting on the boundary at z = Zj is decomposed as (jS.lh where the first term on the rhs is 
the force for a single boundary (when the second one is absent) while the second term is induced 
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by the presence of the other boundary. The interaction part of the force is attractive and it can be 
expressed in the form (|5.2p or, equivalently, as in (j5.4p . The forces acting on the left and on the 
right boundaries are different from each other. For smah separations of the boundaries, as compared 
with the AdS curvature radius, to leading order we recover the result for the Minkowski bulk with 
boundaries. For large separations, the asymptotic expressions for the force are given by (15. 6p and the 
interaction force is again exponentially suppressed. We have also checked with care the transition to 
the Minkowskian limit corresponding to a — t- cxd. For the evaluation of the total vacuum energy in 
the region between the boundaries we have employed zeta function techniques. As the corresponding 
scheme is well described in the literature on the Casimir effect, we have here sketched the main steps 
only. For the total vacuum energy, including the contributions coming from the regions z < zi and 
z > Z2-, one has (j5.16p with the interference part being given by (j5.17p . In the second part of Sect. [5] 
we have described in detail how from our results the corresponding formulae for fermionic Casimir 
densities in Randall-Sundrum-type braneworld scenarios immediately follow. 
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